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spinor polaritons, e.g., Refs. I13U14L in the following, we 
generalize this approach to the case of multiple species of 
excitons, i.e., paraexcitons and orthoexcitons, adopting 
a mean-field coupling scheme between the components^ 
The multi-component exciton gas is considered in sec- 
ond quantization. We start from the Hamiltonian for a 
if-component system in the grand canonical ensemble: 
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The system under consideration is a multi-component gas of interacting paraexcitons and orthoex- 
citons confined in a three dimensional potential trap. We calculate the spatially resolved optical 
emission spectrum due to interband transitions involving weak direct and phonon mediated exciton- 
photon interactions. For each component, the occurrence of a Bose-Einstein condensate changes the 
spectrum in a characteristic way so that it directly reflects the constant chemical potential of the 
excitons and the renormalization of the quasiparticle excitation spectrum. Moreover, the interaction 
between the components leads, in dependence on temperature and particle number, to modifications 
of the spectra indicating phase separation of the subsystems. Typical examples of density profiles 
and luminescence spectra of ground-state paraexcitons and orthoexcitons in CU2O are given. 

PACS numbers: 78.20.-e,78.30.-j,71.35.Lk 



I. INTRODUCTION 

Excitons in semiconductors have been promising candi- 
dates for the observation of Bose-Einstein condensation 
for several decades. At present, cuprous oxide (CU2O) 
is in the focus of experimental efforts due to the large 
binding energy and long lifetime of the exciton states. In 
order to obtain sufficiently high densities, entrapment by 
an external potential is an approved method. 

The theoretical description of excitons in potential 
traps has been carried out so far mostly in the frame 
of a model of ideal bosons. Concepts for the inclu- 
sion of the interaction are well known from the theory 
of atomic condensates ji— and first applications to exci- 
tons exist, too.— Recent investigations in the framework 
of a mean-field formalism in local density approxima- 
tion have shown distinct signatures of a condensate in 
the decay luminescence spectrum of the non-condensed 
excitons.— It is the aim of the present paper to intro- 
duce a generalization of this theory to a multi-component 
gas of interacting paraexcitons and orthoexcitons, where 
the consequences of the interaction on the condensation 
process are of particular interest. We show results for 
the densities of the individual components and their spa- 
tially resolved luminescence spectra for several parameter 
regimes and highlight experimentally relevant cases. 



II. THERMODYNAMICS OF EXCITONS IN A 
POTENTIAL TRAP 

The thermodynamics of a one-component Bose gas has 
been investigated in detail, see, e.g., Refs. LU-|4]. First ap- 
plications of these concepts to excitons have been pre- 
sented in Ref. Jjj and, looking at spectral signatures of 
a condensate, in Ref. |Q. In analogy to generalizations 
for multi-component atomic gases, e.g., Refs. I7H13 and 
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with respective external potentials Vi and chemical po- 
tentials [ii for each species. We assume a contact po- 
tential for the exciton-exciton interaction, with the ma- 
trix hij containing the intra- and inter-species interaction 
strengths. Its components are given by the s-wave scat- 
tering lengths a*-: 
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The Bose field operator ipi obeys the Heisenberg equation 
of motion 
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We decompose the field operators ipi in the usual fashion, 

V>i(r,t)=* 4 (r) + &(r,t), (4) 

where $>i is the (scalar) condensate wave function with 
$i(r) — (i/)i(r, t)) = (ipi(r)) and V^i is the operator of 



the thermal excitons. Inserting the decomposition © 
into Eq. © and following the steps of Ref. [H we obtain 
(arguments dropped for the sake of brevity) 
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with ny = <?*^i +ftjj, wiy = $j<Pj + niy, and the normal 

and anomalous averages riy = (iplipj) and my = (ipiifrj), 
respectively. Equation ([5]) generalizes the familiar Gross- 
Pitaevskii equation (GPE) by including (i) the coupling 
to the thermal excitons and (ii) the coupling of multiple 
components. 

In a first approximation, we neglect all non-diagonal 
averages, i.e., my = ™ij = m y = n ij =0 Vi 7^ j, 
reducing ([SJ and © to effective one-species equations 
with an additional mean field contribution from the other 
species: 
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Thanks to this simplification, Eq. © can be formally 
solved by a Bogoliubov transformation, 

£, = £ h'W^We-*'* + «; (a)aj (o-Je^W'/' 

(9) 
where a enumerates the quasiparticle states. Thereby 
the Bogoliubov amplitudes u, and Vi satisfy the relation 
J2a[ u i(< 7 ) 2 ~ v i( a ) 2 ) — 1- The excitation spectrum -Ej(cr) 
is given by the solution of the eigenvalue problem 
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Equations © and (fTOj) are the multi-component gener- 
alizations of the Hartree-Fock-Bogoliubov (HFB) equa- 
tions. They represent a system of 3K equations which 
are coupled via the interaction matrix elements /jy ■ 

To guarantee gapless spectra, we next apply the Popov 
approximation, i.e., we neglect the anomalous averages 
fhu in © and (ITU|) and find 
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Since the extension of the potential trap is large com- 
pared to the typical length scale of the system (e.g., the 
thermal deBroglie wavelength of the excitons), we can use 
the local density approximation (LDA). Then the exci- 
tons are treated as a locally homogeneous system and 
the spatial dependence enters only via the trap poten- 
tial. In that case, the Bogoliubov equations are readily 
solved, yielding the density nf = nu of thermally excited 
excitons as 
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with ub(E) = [exp(E/kBT) — 1] 1 being the usual Bosc 
function. The excitation spectrum Ei is explicitly given 
by 
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with n\ = |^i| 2 and m = nu — n J + n t- In consistence 
with the LDA, we apply the Thomas-Fermi approxima- 
tion to the Gross-Pitaevskii equation, neglecting the ki- 
netic energy term in fj 12[) . Then we obtain finally for the 
densities of the condensates 
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The expressions P^)) - pT)) have to be solved self- 
consistently Although they look very similar to the one- 
component case, a coupling between the components ap- 
pears via L t and n?. 

In what follows we calculate the densities of excitons 
in CU2O in a strain induced potential trap. 16 In addition 



to the paraexcitons (labeled 'p'), two spin projections of 
orthoexcitons are captured by the trap, denoted by '+' 
and '— ', while the zero component is expelled and plays 
no role. Thus, the number of components K — 3. In 
addition to the usual symmetry of the interaction ma- 
trix, hij — hji with 
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-, in our case it holds 
leaving four indepen- 
As exten- 
sive works on two-component systems^— ^ have shown, 
one of the most interesting aspects of multicomponent 
systems — the occurrence of phase separation — is closely 
tied to the proportions of inter- and intra-species inter- 
action strengths. 

According to Eq. ©, the interaction strengths are 
given by the s-wave scattering lengths of the correspond- 
ing channels, which can be obtained by the solution 
of the four-particle scattering problem. The case of 
positronium-positronium scattering some time ago re- 
ceived much attention 1 ^— and quite reliable values of 
the scattering length for both the singlet and the triplet 
channel have been obtained. In contrast, the description 
of exciton-exciton interaction is a long-standing prob- 
lem, and so far no satisfying solution for the general 
case has been obtained. Especially for CU2O, we ex- 
pect a strong effect of the non-parabolicity of the valence 
band 2 - and of the rather large electron-hole exchange 
interaction 2 ^ on the scattering lengths. Therefore, the 
values we deduced from the scattering lengths of the 
positronium problem given in Ref. [l7| (h++ =0.71 h pp , 

h p+ = 0.33 h pp and h^ = 1.77 h pp with h pp = 7.5 x 10~ 4 

ueV um 3 ), should be considered as representative, only. 
Nevertheless, they allow to show the general behaviour 
of the multi-component exciton system. In the calcu- 
lation we also neglect the difference in the paraexciton 
and orthoexciton mass due to the fc-dependent exchange 
interaction. 22 



III. LUMINESCENCE SPECTRUM 

Excitons decay by emitting photons. This takes place 
either directly, whereby momentum conservation requires 
that only excitons with the same momentum as the emit- 
ted photons are involved, or with assistance of momen- 
tum supplying phonons such that all exciton states can 
participate in the optical emission. Because the opti- 
cal wavelength of the emission is much smaller than the 
trapped exciton cloud, we apply a local approximation 
to the emission spectrum, which, for the homogeneous 
case, is determined by the excitonic spectral function 
A(k,w)2i24: 

Ji(r,«) oc 2^(k = 0)| 2 <5(fiw'-/i,K(r) (18) 

+ y^ \Si(k)\ 2 n B (h\u' - /j, i )A i (r,'k,fkv' - Hi) 

with Si(k) representing the exciton-photon coupling. 
The spectral function is given by the Bogoliubov ampli- 
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In order to account for a finite spectral resolution be- 
ing important for comparison with measured spectra, we 
convolute the spectral intensity (fT5|) with a slit function 
of the shape exp[— (w/A) 4 ]. Here, A is a measure for 
the spectral resolution. Furthermore, in a typical experi- 
mental situation, one images a small stripe of width 2Ax 
elongated along the z-direction onto the entrance slit of 
a spectrograph. Integrating over the y-direction perpen- 
dicular to z we obtain the spatially resolved spectrum of 
the thermal excitons (in the following, the direct conden- 
sate contribution is not considered): 
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with e± (u/, k, r) = (flu' - /j, ± E(k, r))/A. 

In the case of phonon-assisted transitions (being rele- 
vant for the orthoexcitons), we have u' = w — E g x/h — 
Wphonon with E g x being the excitonic band gap of the 
semiconductor. We assume 5(k) to be a constant. Then 
the first term in Eq. (1181) gives rise to a 5-shaped lumi- 
nescence line at the position of the chemical potential, 
with a strength determined by the condensate density. 

For trapped paraexcitons, the zero-phonon decay is rel- 
evant and can be treated by setting u/ = u> — E g x /h and 
S(k) — So8(k — ko). Here ko is the wavevector of the 
intersection of photon and exciton dispersion. Its modu- 
lus is given by fco — E g xn/hc, where n is the refraction 
index and c is the vacuum velocity of light. Due to the 
form of 5(k), the condensate itself does not contribute 
to the direct luminescence process. 

However, as discussed for the one-component exciton 
gas in Ref. @, in case of a potential trap there will be 
indirect signatures of the condensate in the spatially re- 
solved luminescence spectrum. The spectral lineshape 
follows the density distribution in the trap, which in turn 
is bordered by the minimal excitation energy E(k = 0, r). 
For a noncondensed gas the latter quantity is roughly 
parabolic due to the trapping potential, while it is zero in 
the presence of a condensate. Thus with increasing par- 
ticle number (or decreasing temperature) the flat bottom 
of the spectrum of thermal excitons may be a footprint 
of Bose-Einstein condensation. 



IV. RESULTS 

We evaluate the density distributions of the trapped 
excitons in an iterative way. In each step, we keep the 
distributions of two of the components fixed. Under the 
constraint of fixed particle number, we iterate the subset 
of equations ([T3])-([T7]) belonging to the third component 
to self-consistency. We cycle through the components 
until self-consistency of the whole system (ri4"]) - ([TT]) is 
achieved. 
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FIG. 1: Potentials and density profiles in ^-direction at 
(x, y) — (0, 0), and luminescence spectra for paraexcitons (left 
column) and ortho( — )excitons (right column) for a tempera- 
ture of T = 2.2 K and particle numbers of N p = 5 x 10 9 and 
N- — N+ — 5 x 10 in the trap. The corresponding chemical 
potentials are (m p = —2260 |J.eV and fj,- = /j, + — —5920 [ieV. 
Upper row: external trap potential V%, quasiparticle en- 
ergy at k = shifted by fi (i.e., renormalized potential) 
E(0,z,0), and the same quantity without interspecies inter- 
action Eid(0, z, 0). Middle row: densities of thermal excitons 
with (n T ) and without interspecies interaction (nf d ). Lower 
row: luminescence spectra. 



paraexcitons, (iii) only one species of orthoexcitons, (iv) 
both species of orthoexcitons, (v) paraexcitons and one 
species of orthoexcitons, and (vi) all species. To analyze 
these cases, we set the particle numbers of each compo- 
nent to one of two values: N{ — 5 x 10 9 or Ni = 5 x 10 8 for 
i = p, +, — , respectively. We get a rough estimate of the 
corresponding critical temperatures by applying a har- 
monic approximation to the Hertzian potentials. Then, 
a simple Thomas- Fermi calculation for the single compo- 
nent case2£ yields T c ° w 2 K for N = 5 x 10 9 and T c ° w 1 K 
for N = 5 x 10 8 . 

For our calculations, we use values of A = 41 ueV for 
the spectral resolution and Aa; = 25 um for the entrance 
slit of the spectrograph being typical for a triple high 
resolution spectrograph used in the current experiments 
which are under way^ 




FIG. 2: Same presentation as in Fig.[T]but for a temperature 
of T — 1.2 K and particle numbers of N p — 5 x 10 9 and 
N- = N+ = 5 x 10 8 in the trap. The chemical potentials 
are /x p = —2140 (i.eV and /i_ = /i+ = —5610 (i.eV. In the 
middle row, additionally the densities of condensed excitons 
with (n c ) and without interspecies interaction (n1 d ) appear. 



Depending on the temperature and particle numbers 
in the trap, there may occur six distinct situations, fea- 
turing a condensate of (i) none of the species, (ii) only the 



In the following figures we show the respective trap po- 
tentials Vi, minimal excitation energies Ei(k. = 0, z, g = 
0), and density distributions of paraexcitons and orthoex- 



citons versus the z-coordinate. We compare the results 
to the case without inter-component interaction (labeled 
'id'). 



0.4 






> °- 2 

CD 

I, 0.0 


"V^; 


l^S 


Q. 

J, -0.2 




— ■V(z,0)-|i 


-0.4 




— E(0,z,0) 
-.E id (0.z.O) 




FIG. 3: Same presentation as in Fig. [T] but for a temperature 
of T = 1.2 K and particle numbers of N p = N+ = 5 x 10 8 , 
and N- = 5 x 10 in the trap. The chemical potentials are 
H P = -2310 ueV, fi+ = -5550 ueV, and p_ = -5490 ^eV. 



We first investigate case (i) by setting N p — 5 x 10 9 , 
N± = 5 x 10 8 and a temperature T = 2.2 K well 
above both of the estimated critical values. As Fig. Q] 
shows, in the absence of any condensate, the lineshapes 
2?j(k = 0, z, g = 0) roughly follow the external poten- 
tials and the presence of multiple components causes only 
a weak additional renormalization. Both para- and or- 
thoexciton densities — the latter being equal for '+' and 
' — ' species — concentrate in the centers of their traps. A 
noticeable redistribution of the orthoexcitons with re- 
spect to the one-component case results from the large 
number of paraexcitons as well as from the ortho-ortho 
interaction. The corresponding luminescence spectra of 
thermal excitons are shown in Fig. [TJ lower row. Be- 
cause the modulus of the photon vector |k | « 30 urn -1 
is rather small, the integrated zero-phonon spectrum of 
the paraexcitons almost directly resembles the minimal 
excitation energy E p (k = 0, z, g — 0). In the case of 



orthoexcitons, every k- vector contributes and we find a 
broad energy distribution above E±(k = 0, z, g = 0). 

Keeping the particle numbers constant, we lower the 
temperature to T = 1.2 K and show case (ii) in Fig. 
[2] Now the renormalized potential of the paraexcitons 
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FIG. 4: Same presentation as in Fig.[T]but for a temperature 
of T = 1.2 K and particle numbers of N p - 5 x 10 8 , N- = 
iV+ = 5 x 10 9 in the trap. The chemical potentials are fj, p = 
-2300 ^eV, ft- = -5420 ^eV, and fi + = -5470 ueV. 

(Fig. [2l upper left panel) is cut at the chemical potential 
causing an almost fiat bottom of the luminescence spec- 
trum (Fig. [2 lower left panel). Again the densities of 
thermal and condensed paraexcitons show no significant 
deviation from the single component case (Fig. [2j mid- 
dle left panel). In contrast, while isolated orthoexcitons 
would have been condensed, there is no ortho-condensate 
in the fully interacting case, which shows that the pres- 
ence of multiple repulsive components lowers the critical 
temperature. Due to the even higher concentration of 
paraexcitons and the different minimum positions of the 
external potentials (z p — 174 jim and z± = 164 pm), the 
thermal orthoexcitons are slightly pushed aside (Fig. [2] 
middle right panel) . Their spectrum (Fig. [5J lower right 
panel) is qualitatively nearly unchanged with respect to 
case (i). However, due to the lower temperature, the 



spectrum is less widespread. Furthermore, the chemical 
potential nearly touches the renormalized potential caus- 
ing already a smoother curvature of the spectral shape. 
If we exchange the particle numbers of paraexcitons 
and one species of orthoexcitons, i.e. N + = 5 x 10 9 , 
N p = N- = 5 x 10 8 , we realize case (hi), which is pre- 
sented in Fig. [3] While in this case the density distribu- 
tions of thermal and condensed ortho(+)excitons deviate 
only weakly from the one-component case, the displace- 
ment of the thermal paraexcitons is expressed in a heavily 
distorted zero-phonon spectrum. 
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FIG. 5: Same presentation as in Fig.[T]but for a temperature 
of T = 1.05 K and particle numbers of N p = 5 x 10 9 , N- = 
N+ = 5 x 10 8 in the trap. The chemical potentials are fi p = 
-2140 peV, n- = -5585 peV, and fi+ = -5580 peV. 

Increasing also the particle number of the remaining 
ortho-species by an order of magnitude, we generate case 
(iv) , depicted in Fig. [4] As Shi et al£ showed, even for 
finite temperature the condition for phase-separation of 
mutually interacting trapped condensates coincides with 

the T = result of Ho and Shenoy>£ h\_ > h ++ h 

Due to their strong repulsion, the two ortho-condensates 
fulfill this condition and separate into a ball-and-shell 
structure with finite overlap, as seen in Refs. and §. 
Yet, as found in Ref. QJ, at T > no pure '+'- or '-'- 



phases exist and the respective thermal particles are not 
entirely expelled. References @ and [27] pointed out that 
in general the component with the weaker self-interaction 
forms the outer shell. In the case of orthoexcitons, this 
criterion does not apply and the labels '+' and ' — ' can 
be interchanged in all the results presented here. In prin- 
ciple, there should exist an unstable solution with equal 
distributions of the ortho-condensates. The breaking of 
this symmetry is a consequence of our iterative numer- 
ical method. Because of the symmetry of the interac- 
tion, the paraexcitons react to the combined density of 
the orthoexciton species. That is why the distortion of 
the para-spectrum is strongest in the area of overlapping 
ortho-condensates. 
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FIG. 6: Upper and lower row: same presentation as in Fig.[T] 
middle row: density profiles and potentials for paraexcitons 
(left column) and both ortho-species (right column) for a tem- 
perature of T = 0.8 K and particle numbers of N p = 5 x 10 9 , 



iV_ = JV+ = 5 x 10 s 
Up = -2140 ueV, fi~ 



in the trap. The chemical potentials are 
= -5560 [ieV, and fi+ = -5570 ueV. 



N p = 



Let us switch back to the parameters of case (ii), i.e., 
5 x 10 9 and N + = 7V_ = 5 x 10 8 , and lower 
the temperature to T = 1.05 K. As Fig. [5] shows, be- 
sides condensed paraexcitons we now find a small con- 
densate of only one species of orthoexcitons (case v). 



While in the single-component case, for N + — N_, both 
ortho-species have the same critical temperature, now 
the mutual repulsion prevents the simultaneous conden- 
sation of the second species. We have to lower the tem- 
perature to T = 0.8 K, to get condensates of all the 
components (case (vi), Fig. HJ). While the ortho-species 
again form a ball-and-shell structure, ortho- and para- 
condensates do not separate because of their weak inter- 
action (hp + < h pp h ++ ). 



V. CONCLUSION AND OUTLOOK 

We have presented a theoretical approach for the de- 
scription of multi-component interacting excitonic gases 
in potential traps. The resulting system of equations 
has been subject to a number of approximations to 
make it numerically feasible. Finally, coupled multi- 
component equations for the densities of thermal ex- 
citons in Hartrec-Fock-Bogoliubov-Popov approxima- 
tion and the condensate densities following from the 
Gross-Pitaevskii equation in Thomas-Fermi approxima- 
tion have been obtained and numerically solved. Com- 
pared to previous calculations^ an experimentally realis- 
tic, anharmonic trap potential has been used. Six "typ- 
ical" (but not necessarily experimentally realizable) sit- 
uations leading to BEC in one or more of the species 
have been compared. The spatially resolved decay lumi- 
nescence spectra of thermal paraexcitons and orthoexci- 
tons exhibit clear signatures of a condensate. On the one 
hand, there is a flat bottom at the chemical potential 
known from the single-component case£ On the other 
hand, the interspecies interaction causes additional mod- 
ifications of the spectra: if there is a condensate in one 
of the species, the spectrum of the respective other com- 
ponent is clearly distorted. 

In a typical experiment in the bulk or involving rather 
shallow potential traps, orthoexcitons are produced by 



laser excitation, but are converted fast into paraexcitons 
at a rate of 0.3 ns -1 ! 28 ' 29 Therefore, under quasicqui- 
librium conditions, the particle number ratio ortho- 
/paraexcitons is small which corresponds to the cases 
(i), (ii), (v), and (vi). By increasing the stress, however, 
the conversion rate decreases by more than an order of 
magnitude^ By continuous excitation of orthoexcitons, 
therefore, it should be possible to obtain higher orthoex- 
citon/paraexciton ratios like in cases (iii) and (iv). For 
an exciton number of 5xl0 9 , the density in the center 
of the trap is about 10 17 cm" 3 which is experimentally 
achievable. The same holds for the considered tempera- 
tures of T = 0.8. ..2.2 K^ 

Tabic [J summarizes the essential information obtained 
from the cases (i)-(vi) discussed above. Obviously, three 
conclusions can be drawn: First, if the temperature is low 
enough (below the respective critical temperature) , every 
species can form a condensate. Its primary signature is a 
flat bottom of the respective spectrum. Second, a spatial 
separation occurs only between the condensates of the 
two ortho species, because of their strong repulsion. It 
shows up only in the densities, not in the combined spec- 
trum. Therefore, it is important for future experiments 
to measure also the spectrally integrated density profile. 
Third, at occurrence of any ortho but no para conden- 
sate, the para spectrum is distorted in a characteristic 
way indicating a condensate in at least one of the other 
species. 

The presented theory is obviously only a first step to- 
wards a deeper understanding of the physics of trapped 
excitons. To go beyond that includes the solution of the 
Gross-Pitaevskii equation (without the Thomas-Fermi 
approximation) and the inclusion of anomalous densities 
already on the single-component level. Moreover, a gen- 
eral multi-component theory requires the consideration 
of mixed averages which overcomes the effective single- 
component picture but complicates the Bogoliubov trans- 
formation remarkably. 
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Case 


Temperature 










Orthoexciton 


Paraspectrum 


No. 


(K) 


iV p /5 x 10 9 


JV+/5 x 10 9 


JV-/5 x 10 9 


Condensate? 


separation? 


deformation? 


0) 


2.2 


1 


0.1 


0.1 


- 


- 


- 


(ii) 


1.2 


1 


0.1 


0.1 


P 


- 


- 


(iii) 


1.2 


0.1 


0.1 


1 


O- 


- 


/ 


(iv) 


1.2 


0.1 


1 


1 


o+,o- 


/ 


/ 


(v) 


1.05 


1 


0.1 


0.1 


p,°- 


- 


- 


(vi) 


0.8 


1 


0.1 


0.1 


p,o+,o- 


/ 


- 



TABLE I: Summary of cases (i)-(vi): occurrence of a condensate in one or more species, of spatial separation between ortho(+) 
(o+) and ortho( — ) (o_) excitons and of a deformation of the thermal para exciton spectrum in dependence on temperature 
and particle numbers 
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